Generation of modal and path entangled photons 
using domain engineered integrated optical waveguide device 
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Integrated optical devices are expected to play a promising role in the field of quantum information 
science and technology. In this paper, we propose a scheme for the generation of non-degenerate, 
co-polarized, modal and path entangled photons using a directional coupler and an asymmetric Y- 
coupler geometry in type phase-matched, domain engineered lithium niobate (LN) waveguide. The 
nonlinearity in LN is tailored in such a way that quasi phase matching conditions for two different 
spontaneous parametric down conversion (SPDC) processes are obeyed simultaneously, leading to 
a modal and path-entangled state at the output. Assuming typical values of various parameters, 
we show, through numerical simulations, that an almost maximally entangled state is achievable 
over a wide range of waveguide parameters. For the degenerate case, the proposed scheme gives a 
NOON state for N = 2. The generated entangled photon pairs should have potential applications in 
quantum information schemes and also in quantum metrology. By appropriate domain engineering 
and component designing, the idea can be further extended to generate a hyperentangled and two- 
photon multi-path entangled states which may have further applications in quantum computation 
protocols. 

PACS numbers: 



I. INTRODUCTION 

Quantum entangled photons, owing to their unique 
properties, play a crucial role in the implementation 
of quantum information protocols, quantum teleporta- 
tion [U-IH, and quantum computing [4|, @. The entan- 
glement of photons can be exhibited in terms of their 
polarization which is intrinsically binary, or spatial and 
spectral degrees of freedom which are continuous vari- 
ables. An important class of entangled photon states is 
NOON states or path-entangled states, defined as:\N) = 

(\N, 0) o b + |0, N) a jjwhich is the coherent superpo- 
sition of N photons in two different spatial modes a and 
b. Nowadays, such NOON states are being exploited to 
beat the classical Rayleigh diffraction limit and thus, be- 
lieved to find potential applications in quantum metrol- 
ogy, imaging and lithographic techniques 

Spontaneous parametric down conversion (SPDC) pro- 
cess in second order (x 2 ) nonlinear bulk crystals and 
waveguides is a well established method for the gener- 
ation of polarization entangled [9l-[l2| and time-bin en- 
tangled photons [13|,llJ|. The use of waveguide structures 
leads to enhanced nonlinear efficiencies due to tight con- 
finement of the interacting waves and in addition, pro- 
vides a means to control the spatial characteristics of 
the down converted photons by confining them to well 
defined discrete transverse spatial modes [ID, Hrj |. This 
is in contrast to the continuous spatial distribution gov- 
erned by momentum conservation in bulk crystals. By 
restricting the modes of the waveguide structure to the 
two lowest order modes, the mode number can serve as a 
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modal qubit [l7T |. which can be used as an alternative to 
polarization qubit. The guided- wave downconverted pho- 
ton pairs can thus be emitted in the two allowed modes 
leading to modal entanglement. Such waveguide struc- 
tures are also compatible with integrated optics technol- 
ogy and form basic components in photonic circuit de- 
signs, which are expected to play an important role in 
quantum information technology, owing to the compact- 
ness of the device and low loss. The generation, separa- 
tion and processing of entangled photons can be carried 
out on a single chip and hence, the entangled photon 
state is less susceptible to decoherence. The photonic 
circuits based on titanium indiffused channel LN waveg- 
uides have also been recently investigated for different 
combinations of modal, spectral and polarization entan- 
gled photon pairs [lg, LL9| • 

In this paper, we address the issue of generation of 
a modal-path entangled state using SPDC in a domain 
engineered integrated optical waveguide device. Sev- 
eral schemes have been proposed over the years to gen- 
erate two photon and multiphoton path entanglement 
using SPDC process with the use of interferometric 
set ups consisting of polarization beam splitters (PBS) 
or Mach Zehnder interferometers (MZI) |2fj| - [22l |. cavity 
QED [|23(. We propose a scheme for the generation of 
non-degenerate, co-polarized modal entangled state us- 
ing a directional coupler in type phase-matched doubly 
periodically polarized lithium niobate (LN) waveguide. 
The nonlinearity in LN is tailored in such a way that 
quasi phase matching (QPM) conditions for two different 
SPDC processes are obeyed simultaneously (24|. These 
two SPDC processes, when enabled together, lead to a 
mode entangled output state. We have shown through 
quantum mechanical analysis and numerical simulations, 
that for our choice of waveguide geometry and the modes, 
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maximally mode entangled photon pairs can be achieved 
for a wide range of input parameters of the waveguide. 
Using an asymmetric Y-coupler at the output, it is then 
possible to convert the mode entangled state to a path 
entangled state. The use of the directional coupler in 
our scheme facilitates the choice of the modes required, 
which leads to an appreciable overlap of the fields at the 
pump, signal and idler wavelengths. The overlap inte- 
grals, which primarily determine the down conversion 
efficiency, are quite high in comparison to a two mode 
waveguide, giving rise to a high photon pair generation 
rate. Also, due to the design of the structure, we show 
that an efficient generation of photon pairs is possible 
using the fundamental spatial mode of the pump as com- 
pared to the earlier work in which the first excited spa- 
tial mode had to be chosen in order to achieve good ef- 
ficiencies [17] ■ The proposed device is flexible and the 
waveguide parameters can be modified to generate a hy- 
perentangled state or a multi-path entangled state. 

The organization of the paper is as follows. In Sec. [TT1 
we describe the theory and the design of the integrated 
optical waveguide geometry considered for the genera- 
tion of mode-entangled state and the conversion of it to 
a path-entanglement. In Sec. IIII1 we give a brief account 
of the basic quantum mechanical analysis for the two 
coupled SPDC processes involved in the generation of 
the desired mode-entangled state. In Sec. IIV1 we present 
the results of numerical simulations and quantify entan- 
glement by calculating the von Neumann entropy of the 
output state and show that a maximally entangled state 
is achievable using the proposed device. We also investi- 
gate the emission bandwidths of the two corresponding 
SPDC processes. We conclude with a brief summary in 
Sec. El 



II. WAVEGUIDE GEOMETRY 

We seek to generate a mode entangled state using spon- 
taneous parametric down conversion process in which the 
signal and idler photons get generated in an entangled 
pair of modes. Once generated, such a mode- entangled 
state can be converted to a path-entangled state, such 
that, if one of the down converted photons comes out of 
one of the two (or more) available paths, the other pho- 
ton (generated with the first photon) would also come out 
of the same path and similarly for the other path(s). In 
order to achieve this, we consider a waveguide geometry 
as shown in Fig. [TJ Region I is a single mode waveguide 
at the pump wavelength, followed by region II consist- 
ing of a symmetric Y-splitter. Region III is a directional 
coupler supporting the fundamental (symmetric) normal 
mode and the first excited (antisymmetric) normal mode 
at signal and idler wavelengths. Using an appropriate 
phase matching in the directional coupler region, it is 
possible to generate down converted photons which are 
entangled in mode number. 

The working of the proposed device is as follows: When 
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FIG. 1: (Color online) Waveguide design for generating modal 
and path-entangled photons. Region I is a single mode waveg- 
uide at pump wavelength, region II is a symmetric Y-splitter, 
region III is a directional coupler consisting of doubly periodic 
grating in lithium niobate (LN) substrate (dots and crosses 
represent the polarizing directions in the grating). Region IV 
is an asymmetric Y coupler. Pump, signal and idler are all 
z-polarized, z axis being the optic axis for the LN substrate. 



the fundamental mode at pump enters the input waveg- 
uide, it passes through the symmetric Y-splitter and ex- 
cites the fundamental symmetric normal mode at pump 
wavelength in the directional coupler (region III). The 
directional coupler region is assumed to have a doubly 
periodic domain reversal grating which simultaneously 
satisfies the quasi phase matching (QPM) conditions [24] 
for the following two different SPDC processes: pump 
photon in the fundamental normal spatial mode down 
converting either a) both in the fundamental normal spa- 
tial mode or b) both in the first excited normal spatial 
mode. The substrate used in this region is assumed to be 
a domain engineered LN substrate which results in type 
phase matched down conversion. 

We show in Sec. lIII| that when the QPM conditions for 
both the SPDC processes are satisfied simultaneously, it 
will lead to the generation of a mode entangled signal- 
idler photon pair in the doubly reversed region. Thus, 
in contrast to the standard phase matching based on 
a single process either in the bulk or in the waveguide 
configuration, domain engineering can lead to the direct 
production of co-polarized mode-entangled state without 
the need for an interferometric set up or a 50:50 beam 
splitter . 

To convert the mode-entangled state to a path- 
entangled state, at the output of the domain reversed re- 
gion, we have assumed an asymmetric Y-coupler in region 
IV with both the output waveguides being singlemoded 
and the propagation constant of the upper waveguide be- 
ing higher than the lower waveguide for signal and idler 
frequencies. Thus the fundamental normal modes of sig- 
nal and idler exit from the upper waveguide while the 
first excited normal modes exit from the lower waveg- 
uide leading to a path-entangled state. We show through 
numerical simulations that with an appropriate choice of 
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parameters of the coupler, it is possible to optimize the 
generation of the entangled state. In Sec. IIII1 we present 
a brief account of the quantum mechanical analysis lead- 
ing to the generation of two-photon two-mode entangled 
state in the direction coupler region (region III of Fig. 

ID. 



III. QUANTUM MECHANICAL ANALYSIS 

We consider the process of parametric down conversion 
in a z-cut, x-propagating, domain engineered lithium nio- 
bate waveguide. We assume the pump, signal and idler 
wavelengths to be all having extraordinary polarization 
to make use of the largest non- linear coefficient (CZ33) of 
lithium niobate. For the considered pump powers, the 
pump field is assumed to be a classical field and the signal 
and idler fields are represented by quantum mechanical 
operators. The electric field distributions for pump [fun- 
damental normal mode(O)], signal[fundamental(0) and 
first excited(l) normal modes] and idler[fundamental(0) 
and first excited(l) normal modes) are given by 
Pump(O): 

E P o = ^e pQ (r) E p0 ( e iV>#>*-<->>*) + e -*(ft*>»-<"».*))f (l) 



dependence of the first excited modal fields for signal(s) 
and idler (i) wavelengths, L represents the interaction 
length, e s and e, correspond to the optical dielectric per- 
mittivity of LN for signal and idler; d s o (d s i) and dio 
(fiii) represent the quantum mechanical annihilation op- 
erators corresponding to the fundamental mode (first ex- 
cited mode) of signal and idler. f5 q o (q = p, s, i) corre- 
sponds to the propagation constant for the fundamental 
mode and /3 r i (r = s, i) corresponds to the propaga- 
tion constant for the first excited mode at the respective 
wavelengths. Energy is conserved for both the processes 
such that u) p = lo s -\- uii\ u) p , oj s and u>i being the pump, 
signal and idler frequencies, respectively. 

In order that the output state is entangled, the non- 
linearity of LN substrate in the directional coupler re- 
gion is engineered to satisfy the following two quasi phase 
matching (QPM) conditions, simultaneously, 
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Signal(O): 



where Ai and A2 are the QPM grating periods for the 
individual SPDC processes. The second order non linear 
polarization generated in the medium is given by 



E sQ = 1 J dw s e s0 (r) ^J^^oe^ - al e~^ x )z 

(2) 

Signal(i): 



E sl = i J dLo s e sl (f) J^ia^e^ - a^e"^)* 

(3) 

Idler(O): 



Eio=iJ du^o^^^-ia^e^-ale-^z (4) 
Idler(l): 



E a = i I duiien (r) 



ftlJ. 
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^(aae^-ale-^z (5) 



where the subscripts and 1 refer to the fundamental and 
the first excited modes respectively, e p o(r), e s o (r) and 
e-io (f) represent the transverse dependence of the funda- 
mental modal fields of pump (p), signal(s) and idler(i), 
respectively and e s \ (f) and en (f) represent transverse 
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The indices i, j, k (1, 2, 3) corresponds to x, y, z com- 
ponents of the co-ordinate system. Thus, Ej represents 
the j th component of the total electric field within the 
medium and is the corresponding non-linear coeffi- 
cient. For the case under consideration, the components 
of the total electric field are 



=0, 
E 2 = , 

E3 = EpQ + E s o + E s i + E i0 



(9) 
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The interaction Hamiltonian is derived as [24 
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^33 (Epo E s0 E i0 + E p0 E s iEa) dxdydz 

(10) 

Using the expressions for the electric fields defined in 
Eq. HJni) and under rotating wave approximation (RWA) 
and energy conservation, we obtain the following expres- 
sion for the interaction Hamiltonian: 



4 



Hi n t — / duj s 



\n s0 n i0 J V I 
\n sl n a J V J 



dx(ll) 



Io and I\ denote the overlap integrals of the electric field 
profiles of the pump, signal and idler modes involved in 
the two SPDC processes and are given as 



Io = 



h 



3 P o (r) e s0 (r) e l0 (r) dydz, 



(12) 



e P o (r) e s i (f) en (f) dydz (13) 



The domain engineering in LN substrate results in two 
independent spatial frequency components in the non lin- 
ear coefficient variation, along the propagation direction 
[2^ | . Thus, the effective non linear coefficient d including 
the effect of periodic domain reversal is given as 

t 4c?33 . ij( lX —iK\x iKzx —iK%x\ 



l 33 



-terms at other spatial frequencies (14) 



Replacing ^33 by ^33 in Eq. (jTTj) and simplifying, the 
interaction Hamiltonian in the interaction picture is ob- 
tained as 



H int = I dw s (C , o(aJ d] + a s0 aio) + C 1 (a^ 1 a} 1 + a s ia i i)) 

(15) 
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C n = - I 1 — e 2 smc (16) 
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/4d 33 ^ % /a; s w. i /i£ , po\ . f AkiL\ 



7rn s inii / \ . 

and Afco and Aki are phase mismatch for the two pro- 
cesses given by 

Afc =K X - (3 p0 + Pso + fto (18) 
Aki =K 2 - p p0 + Psi + Pn (19) 
The output two photon state is found as 

du s i(C \so,io} + Ci\s u h)) (20) 

and is entangled; Cq — ~j^C and C± = —j-C x . The 
relative magnitudes of Co and C\ coefficients will deter- 
mine if the state is maximally entangled or not. These, 
in turn, depend upon the overlap integrals and the ef- 
fective indices of the interacting modes at pump, signal 
and idler wavelengths. In Sec. IIV1 we carry out simula- 
tions for a planar waveguide and using practical values 
of waveguide parameters, we investigate various aspects 
of the output entangled state. 



IV. NUMERICAL SIMULATIONS 
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FIG. 2: Waveguide directional coupler geometry correspond- 
ing to Region III of Fig. [T] 



In this section, we report on the results of simulations 
on the entangled states generated in the directional cou- 
pler region using Titanium indiffused domain engineered 
planar lithium niobate waveguide. We model region III as 
shown in Fig. O In order to demonstrate the generation 
of the entangled state in this region, we assume the planar 
waveguides to have step refractive index profiles. Since 
actual devices will use channel waveguides, our analysis 
can be used by first representing the channel waveguides 
by equivalent planar waveguides using the standard ef- 
fective index method pa |. The analysis presented can 
be easily extended to channel waveguide geometries with 
graded refractive index profiles; however, the general con- 
clusions of our analysis would still be valid. As shown in 
Fig- HI the guiding regions (d/2 < \y\ < a) are assumed 
to have refractive index m and the surrounding medium 
to have a refractive index n 2 {< Tlx). Thus, the geometry 
acts like a directional coupler with the cores of the indi- 
vidual waveguides, separated by distance d. The individ- 
ual core widths are (a — d/2) = r. Modal analysis for the 
coupler region was carried out [26] and the eigenvalue 
equations for the symmetric and antisymmetric modes 
were obtained. The solution of these eigenvalue equa- 
tions gives the propagation constants of all the modes, 
at all the wavelengths. For the numerical simulations, 
the refractive indices of the substrate at the correspond- 
ing wavelengths have been calculated u sing temperature 
dependent Sellmeier equation for LN [23|. The pump, 
signal and idler wavelengths are chosen to be 750 nm, 
1452 nm and 1550 nm, the corresponding refractive index 
differences (ni— 112) are assumed to be 0.0033, 0.0026 and 
0.0025, respectively [24,28]. The propagation constants 
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FIG. 3: (Color online) Variation of propagation constant with 
the core separation d corresponding to (a) signal and (b) idler. 



evaluated for the symmetric and antisymmetric modes at 
signal and idler wavelengths vary with the core separa- 
tion d as shown in Fig. [3] We can see from Fig.|3]that, as 
expected, beyond a certain core separation in the coupler 
region, the symmetric and antisymmetric modes have al- 
most the same propagation constants for both signal and 
idler wavelengths. 

The transverse electric field distributions of the inter- 
acting modes at pump, signal and idler are shown in 
Fig. [4] for a core separation of 6 \im. 



Fig. 



4 a 



shows 



the electric field distributions corresponding to the funda- 
mental symmetric modes in the directional coupler region 
at the pump (solid), signal (dashed) and idler (dotted) 
wavelengths and Fig. |4(b)| shows the field distributions 
corresponding to the fundamental symmetric mode of the 
pump (solid) and the first excited antisymmetric modes 
at signal (dashed) and idler (dotted) wavelengths. The 
pump field is more tightly confined than the signal and 
idler wavelengths due to its shorter wavelength. Due to 
the choice of the waveguide geometry and wavelengths, 
it can be seen from the figure that there is a good over- 
lap between the fields of the interacting modes at the 
respective wavelengths, resulting in high down conver- 
sion efficiencies for both the SPDC processes. From the 
field expressions for the spatial modes at pump, signal 
and idler, we have evaluated analytical expressions for 
the overlap integrals Iq and 1\ (defined by Eq. (fT2]) and 
Eq. (|13p ) which are plotted with respect to the core sepa- 
ration in Fig. [5] The figure shows that by an appropriate 
choice of the waveguide separation, it is possible to obtain 
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FIG. 4: (Color online) Normalized transverse electric field 
distributions in the direction coupler region corresponding to 
(a) fundamental symmetric normal modes at pump, signal 
and idler wavelengths, and(b) fundamental symmetric normal 
mode at the pump and first excited antisymmetric normal 
mode at signal and idler wavelengths. 
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FIG. 5: (Color online) Variation of the overlap integrals with 
the core separation d. 



the desired ratios of the overlap integrals that lead to a 
high degree of entanglement between the downconverted 
photons. 

In order to quantify the entanglement, we calculate 
Von Neumann entropy [2£| for the case of perfectly phase 
matched output entangled state of Eq. (PZIH) (after nor- 
malization) as S = —tr(p s log 2 p s ), (tr denoting the trace 
and p s is the reduced density matrix of the signal (s) pho- 
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ton) and is derived as 

S = -J^Jl lo& jjh j^jg log 2 (21) 
The variation of 5" with the waveguide separation is plot- 
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FIG. 6: (Color online) Variation of Von Neumann entropy 
with the core separation d. 

ted in Fig.[SJ Von Neumann entropy (S) for a maximally 
entangled state is 1 and for a product state. We can 
see from Fig. [51 that for the entangled state generated 
using the proposed waveguide device, S becomes almost 
unity beyond a waveguide separation of about 6 /im, thus 
the device is capable of producing a maximally entan- 
gled state over a wide range of waveguide separation. 
The corresponding grating periods required for the two 
SPDC processes at this core separation are 17.521 /zro 
and 17.358 /im respectively, as calculated from Eq. © 
and Eq. (JT]). Another condition for obtaining a maxi- 
mally entangled state with high efficiency requires identi- 
cal bandwidths for the two enabled SPDC processes [24| ■ 
Fig. [7] shows the output spectra at signal wavelength cor- 
responding to both the SPDC processes. It can be seen 
from the figure that the two processes have almost identi- 
cal bandwidths (16 nm) ensuring a maximally entangled 
state over the region of overlap [17j . 
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FIG. 7: (Color online) Normalized output signal spectra cor- 
responding to the two SPDC processes considered. 

We may mention here that the above simulations have 
been carried out using a titanium indiffused waveguide; 




FIG. 8: Waveguide geometry for generating four path entan- 
gled state. 

simulations can also be performed for proton exchange 
(PE) LN waveguides which may further increase the ef- 
ficiencies due to tighter confinement of modes. Thus in 
region III, we achieve two photon modal entanglement. 
In the proposed device, this mode entangled state fur- 
ther gets converted to a path entangled state, when the 
photons reach region IV, consisting of an asymmetric Y- 
coupler where owing to the difference in the propagation 
constants of the upper and the lower waveguides, the fun- 
damental modes of signal and idler would exit from the 
upper waveguide while the first excited modes would exit 
from the lower waveguide. This leads to a path entangled 
state of the form 

\ip)= J du s i(C \s u ,i u ) +d\si,ii) (22) 

where the subscripts u and / corresponds to the upper 
and lower waveguides, respectively. Thus both the sig- 
nal and idler photons would exit either from the upper 
waveguide (fundamental mode) or from the lower waveg- 
uide (first excited mode), resulting in path entanglement. 
If the SPDC process is chosen to lead to a degenerate pair 
of photons, then the output from the device would be a 
NOON state with N=2. By modifying the waveguide 
parameters, the above scheme can be extended to gen- 
erate different combinations of modal, polarization, path 
entangled state giving rise to a hyper-entangled state. 
Also, it can be used to generate a multipath entangled 
state by appropriately engineering the phase matching 
grating in the coupler region and enabling multiple SPDC 
processes simultaneously. This would lead to the gener- 
ation of higher order modes giving rise to a multi-mode, 
two photon entangled state. The output Y coupler can 
then be fabricated with more than two output ports, cor- 
respondingly giving rise to a multi-path, two photon en- 
tangled state. 

The multipath entangled state can also be achieved 
by combining more of these above geometries and de- 
signing the structure appropriately. As an example, a 
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four path-two photon entangled state can be generated 
by designing the structure as shown in Fig. |SJ Two of 
the waveguides of Fig. [1] are combined using a symmetric 
Y-coupler at the input and the signal-idler photon pair 
can exit from either of the four available paths with equal 
probability, giving rise to a maximally multi(four, in this 
case)-path entangled state. The potential applications 
for such multipath entangled states have already been 
discussed in Ref. [30j . 

V. CONCLUSION 

In this paper, we have addressed the issue of genera- 
tion of mode entangled photon pairs using optical waveg- 
uide geometry which has the inherent advantage of being 
compact (making the photons less susceptible to deco- 
herence) and having enhanced optical non linearities in 
comparison to bulk crystal or optical fibers. The pro- 
posed scheme is very efficient as the pump, signal and 
idler are taken to be all extraordinarily polarized, mak- 
ing use of the highest non linear coefficient ^33 of lithium 
niobate, and due to our choice of the geometry and the 



modes, the overlap integrals determining the efficiency 
of the two processes are also significant. We have shown 
that a maximally entangled state is achievable for a wide 
range of core separation of the coupler, with the addi- 
tional advantage that one of the entangled photons is 
emitted at the typical wavelength of 1550 nm and can 
be used as a flying qubit. This mode entangled state 
can further be converted to a path entangled state which 
should have varied applications in quantum lithography 
and quantum information processing. The proposal when 
extended to degenerate SPDC processes would lead to 
the production of a NOON state with N = 2 photons. 
By appropriately choosing the waveguide parameters and 
also designing the integrated structure, the idea is appli- 
cable to produce a hyperentangled state and multi-path 
two photon entangled state which may have future ap- 
plications in quantum computation protocols. 
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